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Abstract 

Jim Propp's rotor router model is a deterministic analogue of a random walk on 
a graph. Instead of distributing chips randomly, each vertex serves its neighbors 
in a fixed order. We analyze the difference between Propp machine and ran- 
dom walk on the infinite two-dimensional grid. It is known that, apart from a 
technicality, independent of the starting configuration, at each time, the number 
of chips on each vertex in the Propp model deviates from the expected number 
of chips in the random walk model by at most a constant. We show that this 
constant is approximately 7.8, if all vertices serve their neighbors in clockwise or 
counterclockwise order and 7.3 otherwise. This result in particular shows that 
the order in which the neighbors are served makes a difference. Our analysis also 
reveals a number of further unexpected properties of the two-dimensional Propp 
machine. 



1 Introduction 

The rotor-router model is a simple deterministic process suggested by Jim Propp. It can 
be viewed as an attempt to derandomize random walks on graphs. So far, the "Propp 
machine" has mainly been regarded on infinite grids Z d . There, each vertex x G Z rf is 
equipped with a "rotor" together with a cyclic permutation (called "rotor sequence" ) of 
the 2d cardinal directions of Z d . While a chip (particle, coin, . . . ) performing a random 
walk leaves a vertex in a random direction, in the Propp model it always goes into the 
direction the rotor is pointing. After a chip is sent, the rotor is rotated according to 
the fixed rotor sequence. This shall ensure that the chips are distributed highly evenly 
among the neighbors. 

The Propp machine has attracted considerable attention recently. It has been shown 
that it closely resembles a random walk in several respects. The first result is due 
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to Levine and Peres [8|, |9j who compared random walk and Propp machine in an 
aggregating model called Internal Diffusion-Limited Aggregation (IDLA) [3]. There, 
each chip starts at the origin of Z d and walks till it reaches an unoccupied site, which 
it then occupies. In the random walk model it is well known that the shape of the 
occupied locations converges to a Euclidean ball in M. d jsj. Recently, Levine and Peres 
@, 0] proved an analogous result for the Propp machine. Surprisingly, the convergence 
seems to be much faster. Kleber 0] showed experimentally that for circular rotor 
sequences after three million chips the radius of the inscribed and circumscribed circle 
differs by approximately 1.61. Hence, the occupied locations almost form a perfect 
circle. Some more results on this aggregating model in two dimensions can be found 
in Section [HI 

Cooper and Spencer [l| compared the Propp machine and the random walk in terms 
of the single vertex discrepancy. Apart from a technicality which we defer to Section [2, 
they place arbitrary numbers of chips on the vertices. Then they run the Propp machine 
on this initial configuration for a certain number of rounds. A round consists of each 
chip (in arbitrary order) doing one move as directed by the Propp machine. For the 
resulting position, for each vertex they compare the number of chips that end up there 
with the expected number of chips that a random walk in same number of rounds would 
have gotten there starting from the initial configuration. Cooper and Spencer showed 
that for all grids Z d , these differences can be bounded by a constant q independent of 
the initial set-up (in particular, the total number of vertices) and the run-time. 

For the case d — 1, that is, the graph being the infinite path, Cooper, Doerr, Spencer, 
and Tardos 0] showed among other results that this constant c\ is approximately 2.29. 
They further proved that to have the discrepancy on a particular vertex maximal it 
suffices that each location has an odd number of chips at at most one time. 

In this paper, we rigorously analyze the Propp machine on the two-dimensional grid Z 2 . 
A particular difference to the one- dimensional case is that now there are two non- 
isomorphic orders in which the four neighbors can be served. The first are clockwise 
and counterclockwise orders of the four cardinal directions. These are called circular 
rotor sequences. All other orders turn the rotor by 180° at one time and are called non- 
circular rotor sequences. We prove C2 ~ 7.83 for circular rotor sequences and C2 ~ 7.29 
otherwise. To the best of our knowledge, this is the first paper showing that the rotor 
sequence can make a difference. 

We also characterize the respective worst-case configurations. In particular, we prove 
that the maximal single vertex discrepancy can only be reached if there are vertices 
which send a number of chips not divisible by four at at least three different times. 
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The remainder of this paper is organized as follows. The basic notations are given in 
Section [21 In Section [H] we roughly speaking show that, by suitably choosing the initial 
configuration, we may prescribe the number of chips on each vertex at each time mod- 
ulo 4. This will yield sharp lower bounds, since in Section H] we see that the discrepancy 
on a vertex can be expressed by exactly this information. In Sections [5] and El we derive 
sufficient information about initial configurations leading to maximal discrepancies on 
a vertex so that we then can estimate the maximum possible discrepancy numerically. 
This estimate is shown to be relatively tight in Section [71 Since the investigation up to 
this point in particular showed that different rotor sequences lead to different results, 
we briefly examine the aggregating model in this respect in Section [HI We summarize 
our results in the last section. 



2 Preliminaries 

To bound the single vertex discrepancy between Propp machine and random walk on 
the two-dimensional grid we need several preliminaries, which will be introduced in 
this section. 

First, it will be useful to use a different representation of the two-dimensional grid Z 2 . 

Let dir : = {(+;), c\),{z{), (;})}■ Define a S ra P h G = via v = {© i 

X\ = X2 (mod 2)} and E = {(x, y) G V 2 | x — y G dir}. Clearly, G is isomorphic 
to the standard two-dimensional grid G' = (Z 2 ,£") with E' = {(x, y) G Z 2 | ||x — 
y||i = 1}. Therefore, our results on G immediately translate to G'. The advantage of 
our representation is that now each direction D G dir can be uniquely expressed as 
D = £z(q) + with e x , s y G { — 1,1}. This allows a convenient computation of the 
probability distribution of the random walk on the grid (see equation (TjQ) below). For 
convenience we will also use the symbols { /*, \, \ } to describe the directions in 
the obvious manner. 

In order to avoid discussing all equations in the expected sense and thereby to simplify 
the presentation, one can treat the expectation of the random walk as a linear ma- 
chine [l|. Here, in each time step a pile of k chips is split evenly, with fc/4 chips going 
to each neighbor. The (possibly non-integral) number of chips at vertex x at time t is 
exactly the expected number of chips in the random walk model. 

For x, y G V and t G No, let x ~ t denote that x± = x-i = t (mod 2) and x ~ y denote 
that x\ = X2 = yi = y2 (mod 2). A vertex x is called even or odd if x ~ or x ~ 1, 
respectively. 
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A configuration describes the current "state" of linear or Propp machine. A configura- 
tion of the linear machine assigns to each vertex x £ V its current (possibly fractional) 
number of chips. A configuration of the Propp machine assigns to each vertex x £ V 
its current (integral) number of chips and the current direction of the rotor. A config- 
uration is called even (odd) if all chips lie on even (odd) vertices. 

As pointed out in the introduction, there is one limitation without which neither the 
results of [l|, 0] nor our results hold. Note that since G is a bipartite graph, chips that 
start on even vertices never mix with those starting on odd vertices. It looks like we 
are playing two games at once. However, this is not true, because chips at different 
parity vertices may affect each other through the rotors. We therefore require the 
initial configuration to have chips only on one parity. Without loss of generality, we 
consider only even initial configurations. 

A random walk on G can be described nicely by its probability density. By i^(x, t) we 
denote the probability that a chip from vertex x arrives at the origin after t random 
steps ( "at time t" ) in a simple random walk. Then, 

i/(x,t)=4-( ( J i)/2 )( (4+ 4 )/2 ) (1) 

for x ~ t and Hx^ < t, and H(x.,t) = otherwise. 

We now describe the Propp machine in detail. First, we define a rotor sequence by 
a cyclic permutation next: dir — > dir. That is, after a chip has been sent in direc- 
tion A, the rotor moves such that afterwards it points in direction next(A). Instead 
of using next directly, it will often be more handy to describe a rotor sequence as 
a 4-tuple TZ = (/ , next(/"), next 2 (/'), next 3 (/')). We distinguish between cir- 
cular and non-circular rotor sequences. Circular rotor sequences are either clockwise 
(/*, \, k/, \) or counter-clockwise (/*, \, fc /, \). All other rotor sequences are called 
non-circular. Our main focus is on the classical Propp mach ine in which all vertices 
have the same rotor sequence. In [lj], Cooper and Spencer allow different rotor se- 



quences for each vertex x. Our results also hold in this general setting. However, to 
simplify the presentation we will typically assume that there is only one rotor sequence 
for all vertices x. 

In the following notations, we implicitly fix the rotor sequence as well as the initial 
configuration (that is, chips on vertices and rotor directions at time t — 0). In one 
step of the Propp machine, each chip does exactly one move, that is, it moves in the 
direction the arrow associated with his current position is pointing and updates the 
arrow direction according to the rotor sequence. Note that the particular order in which 
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the chips move within one step is irrelevant (as long as we do not label the chips). By 
this rule, all subsequent configurations are determined by the initial configuration. For 
all x G V and t G No let /(x, t) denote the number of chips on vertex x and ARR(x, t) 
the direction of the rotor associated with x after t steps of the Propp machine. 

To describe the linear machine we use the same fixed initial configuration as for the 
Propp machine. In one step, each vertex x sends a quarter of its (possibly frac- 
tional) number of chips to each neighbor. Let -E^x, t) denote the number of chips 
at vertex x after t steps of the linear machine. This is equal to the expected num- 
ber of chips at vertex x after a random walk of all chips for t steps. Note that 
E(x, t) = \ £ AeDIR E(x + A, t - 1) by definition. 

3 Mod-4-forcing Theorem 

For a deterministic process like the Propp machine, it is obvious that the initial config- 
uration (that is, the location of each chip and the direction of each rotor), determines 
all subsequent configurations. The following theorem shows a partial converse, namely 
that (roughly speaking) we may prescribe the number of chips modulo 4 on all vertices 
at all times and still find an initial configuration leading to such a game. An analogous 
result for the one- dimensional Propp machine has been shown in |2|. 

Theorem 1 (Mod-4-forcing Theorem). For any initial direction of the rotors and 
any ir: V x No — ■> {0,1,2,3} with vr(x, t) = for all x ^ t, there is an initial 
even configuration with /(x, 0), x G V that results in subsequent configurations with 
/(x, t) = 7r(x, t) (mod 4) for all x and t. 

Proof. Let arr(x, 0) describe the initial rotor directions given in the assumption. The 
sought-after configuration can be found iteratively. We start with /(x, 0) := vr(x, 0) 
chips at location x. 

Now assume that our initial (even) configuration is such that for some T G N we have 
/(x, t) = 7r(x, t) (mod 4) for all t < T. We modify this initial configuration by defining 
/'(x, 0) := /(x, 0) + £ X 4 T for even x, while we have /'(x, 0) = for odd x. Here, 
e x G {0, 1, 2, 3} are to be determined such that /'(x, t) = 7r(x, t) (mod 4) for all t < T. 

Observe that a pile of 4 T chips splits evenly T times. Hence for all choices of the £ x 
we still have /'(x, t) = 7r(x, t) (mod 4) for all t <T. At time T, the extra piles of 4 T 
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chips have spread as follows: 

f(x,T) = /(x,T) + ^{ tZ-vMtZ-vX 

y ^0 ^ 2 /V 2 / 

||y-x||oo<T 

Let initially e y := for all y G V. By induction on ||y||i, we change the e y to their 
final value. We keep e y = for all y with ||y||i < 2T. 

Assume that for some 9 G No, the current e y fulfill /'(x, T) = 7r(x, T) (mod 4) for all x 
with ||x||i < 9. We now determine e y for all y with ||y||i = IT + 9 in such a way that 
f (x,T) = 7r(x,T) (mod 4) for all x G V such that ||x||i < 9. 

Fortunately, to achieve /'(x, T) = 7r(x, T) (mod 4) for some x G V such that ||x||i = 0, 
it suffices to change a single e y , y G V, ||y||i = 2T + 9. Without loss of generality, let 
x G V, ||x||i = 6>, and x ~ T such that Xi, x 2 > 0. Let y = y(x) = (xi+T, x 2 +T). Now 
choosing e y G {0,1,2,3} such that e y = 7r(x, T) — /(x, T) (mod 4) yields /'(x, T) = 
/(x, T) + £ y = 7r(x, T) (mod 4) and /'(x, T) = /(x, T) for all other x 6 ^ such that 

ll x lli < ^« 

Hence for each x 6 V such that ||x||x = ^, we find a y(x) and a value for £ y ( x ) such 
that the resulting /'(x, T) are as desired. All other e y with ||y||i = 2T + 9 remain fixed 
to zero. 

This shows that for all 9 G N, we may choose e y , y G V, such that f'(x,t) = 
7r(x, t) (mod 4) for all t < T and all x G V such that ||x||i < 9. By compactness 
principle, this yields the existence of e y , y G V, such that /'(x, t) = vr(x, t) (mod 4) for 
all t < T and x G 7. 

Up to this point, we proved that for all T G M, there is an even initial configuration 
such that /'(x, t) = 7r(x, t) (mod 4) holds for all £ < T and x G 7. Invoking the 
compactness principle again finishes the proof. □ 

4 The Basic Method 

In this section, we lay the foundations for our analysis of the maximal possible single- 
vertex discrepancy. In particular, we will see that we can determine the contribution 
of a vertex to the discrepancy at another one independent from all other vertices. 
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In the following, we re-use several arguments from [l|, |2|. For the moment, in addition 
to the notations given in Section [21 we also use the following mixed notation. By 
E(~x,ti,t2) we denote the (possibly fractional) number of chips at location x after first 
performing t\ steps with the Propp machine and then t<i — t\ steps with the linear 
machine. 

We are interested in bounding the discrepancies |/(x, t) — E(pc, t)\ for all vertices x and 
all times t. Since we aim at bounds independent of the initial configuration, it suffices 
to regard the vertex x = 0. From 

£(0,0, t) = E(0,t), 
E(0,t,t) = f{0,t), 

we obtain 

t-i 

f{0,t) - E{0,t) = ^2(E(0,s + l,t)-E(0,s,t)). 

s=0 

Now E(0, s + l,t)- E(0, s, t) = E xe y E£ S) + next^ 1 (arr(x, s)), t - s - 1) - 
H(x,t — s)) motivates the definition of the influence of a Propp move (compared to a 
random walk move) from vertex x in direction A on the discrepancy of (t time steps 
later) by 

inf(x, A, t) := H(x + A, t - 1) - #(x, t). 

To finally reduce all arrs involved to the initial arrow settings arr(-,0), we define 
Sj(x) := minju > | i < J^"=o /( x ' ^)} ^ or a ^ * ^ ^o- Hence at time Sj(x) the 
location x is occupied by its i-th chip (where, to be consistent with we start 
counting with the 0-th chip). 

Let T be a time at which we regard the discrepancy at 0. Then the above yields 

/(0,T)-£(0,T) = Yl INF(x,NEXT ? (ARR(x,0)),T- Si (x)). (2) 

xeK i>0, 
«i(x)<T 

Since the inner sum of equation (|2j) will occur frequently in the remainder, let us define 
the contribution of a vertex x to be 

CON(x) := INF(x, next j (arr(x, 0)),T - Sj(x)), 

j>0, 
«i(x)<T 
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where we both suppress the initial configuration leading to the Sj(-) as well as the run- 
time T. Occasionally, we will write CON^ to specify the underlying initial configuration. 

The first main result of this section, summarized in the following theorem, is that it 
suffices to examine each vertex x separately. 

Theorem 2. The discrepancy between Propp machine and linear machine after T time 
steps is the sum of the contributions CON(x) of all vertices x, i.e., 

f(0,T) - E(0,T) = ]Tcon(x). 

x€V 

Our aim in this paper is to prove a sharp upper bound for the single- vertex discrepancies 
|/(y, T) — E(y, T)\ for all y and T. As discussed already, by symmetry we may always 
assume x = 0. To get rid of the dependency of T, let us define maxcon(x) to be 
the supremum contribution of x over all initial configurations and all T. We will 
shortly see that the supremum actually is a maximum (Corollary [TOl). that is, there 
is an initial configuration and a time T such that CON(x) = maxcon(x). Since the 
contribution only depends on T — Sj(x) and the (Mod-4)-forcing theorem tells us how to 
manipulate the Sj(x), we may choose T as large as we like (and still have a configuration 
leading to CON(x) = maxcon(x)). Provided that Xlxev maxcon(x) is finite (which 
we prove in the remainder), we obtain that ^ xe y maxcon(x) is a tight upper bound 
for sup(/(0, T) — E(0, T)), where the supremum is taken over all initial configurations 
and all T. 

To bound |/(0, T) — E(0,T)\, we need an analogous discussion for negative contribu- 
tions. Let mincon(x) be the infimum contribution of x over all initial configurations 
and all T. Fortunately, using symmetries, we can show that Yl x &v MAX CON(x) = 
— Sxev mincon(x), hence it suffices to regard positive contributions. Let us shortly 
sketch the symmetry argument and then summarize the above discussion. 

Observe that sending one chip in each direction at the same time does not change 
CON(x). That is, for all x and t we have 

^lNF(x,A,t) = 0. (3) 

Agdir 

This follows right from the definition of INF and the elementary fact H(x,t) = 
I X^AeDiR ^( x + A,t — 1). Based on equation we will ignore piles of four chips 
(and multiples) at a common time t in the remainder of this section. The remaining 
one to three chips are called odd chips. Note that there is no relation between odd 
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chips and odd vertices/configurations as denned in Section [21 

To describe the symmetries of CON, we further distinguish the non-circular ro- 
tor sequences. We call (/\\,\,,/) and (/\^,\,\) x-alternating and (f, \ 
,\,y/) and (/*, / , \, \) y- alternating. Now a short look at the definition of 
maxcon reveals symmetries like maxcon((^)) = maxcon((~^)) for circular ro- 
tor sequences, maxcon((^)) = MAXCON((_f^ )) for x-alternating rotor sequences, and 
MAXCON ((^)) = maxcon((~^ 1 )) for y-alternating rotor sequences. The following 
lemma exhibits symmetries for maxcon and mincon. It shows that the discrepancies 
caused by having too few or too many chips have the same absolute value. 

Lemma 3. For all x G V , the following symmetries hold for 

• circular rotor sequences: MAXCON((^)) = — MINCON((~^ 1 )), 

• x-alternating rotor sequences: MAXCON ((^)) = — MINCON((~^ 1 )), 

• y-alternating rotor sequences: MAXCON = -MINCON ((_?£ ))■ 

Proof. The proofs are not difficult, so we only give the one for the first statement. 
We show that for each configuration C\ there is another configuration C3 and a simple 
permutation 7r of V with cONd(x) = — CONc 3 (-7r(x)) for all implicit run-times T and 
assuming the clockwise rotor sequence 71 := (/*, \, y, \) for both C\ and C 3 . By The- 
orem!]], there is a configuration C 2 which sends, using the rotor sequence (/*, \, \), 
an odd chip from ("i^ 1 ) in direction at time t if and only if C\ sends an odd chip 

from (^) in direction (^) at time t. Note that CON C2 ((~^)) = CON Cl (x). A configu- 
ration C3 which sends for each single chip C2 sends, three chips from the same vertex 
in the same direction at the same time obeys rotor sequence TZ and gives by equa- 
tion © a contribution CON C3 ((~ X1 )) = -CON C2 ((~^ 1 )) = -CON Cl (x). In consequence, 
mincon((~^)) = — maxcon((^)) for the clockwise rotor sequence TZ. □ 

Now Lemma [3] immediately yields ^ xgy mincon(x) = — ^ xgy maxcon(x). There- 
fore, it suffices to regard maximal contributions. 
Theorem 4. 

sup |/(0, T) - E(0,T)\ = V maxcon (x) 
is a tight upper bound for the single vertex discrepancies. 
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5 The Modes of inf 



In Theorem H] we expressed the discrepancy as sum of contributions con(x), which 
in turn are sums of the influences inf(x, A, £). To bound the discrepancy, we are 
now interested in the extremal values of such sums. In this section we derive some 
monotonicity properties of these sums. For this, we define 

INF(x, A,t) := ^2 inf(x, A,t) 
Ae^4 

for a finite sequence A ■— (A^\ A^ 2 \ . . .) of rotor directions ordered according to a fixed 
rotor sequence. In the remainder of the article all finite sequences of rotor directions 
for which we use the calligraphic A are ordered according to their respective rotor 
sequence. 

Let X C R. We call a mapping /: X — > K unimodal, if there is a t\ G X such that 
f\x<ti as well as /|a;>ti are monotone. We call a mapping /: X — > R bimodal, if there 
are ti,t 2 G X such that f\ x <ti> f\ti<x<t 2 , and f\t 2 <x are monotone. We call a mapping 
/ : X — > M. strictly bimodal, if it is bimodal, but not unimodal. In the following, we 
show that all inf(x, A, t) are bimodal in t. 

From equation ([3]) we see that 



lNF(x,(A (1) ,A (2) ,A (3) ),t) = -lNF(x,DlR\{A (1) ,A (2) ,A (3) },t) and ^ 
inf(x, (A (1) ,...,A (fe) ),t) =inf(x, (A (1) ,...,A (fc " 4) ),t) for k > 4. 

This shows that it suffices to examine inf(x, A, t) for A of length one and two, which 
is done in Lemmas [6] and d respectively. For both proofs, we need Descartes' Rule of 



Signs, which can be found in H 



Theorem 5 (Descartes' Rule of Signs) . The number of positive roots counting multi- 
plicities of a non-zero polynomial with real coefficients is either equal to its number of 
coefficient sign variations (i.e., the number of sign changes between consecutive nonzero 
coefficients) or else is less than this number by an even integer. 



With this, we are now well equipped to analyze the monotonicity of INf(x, A, •) for 

|^| e {1,2}. 

Lemma 6. For all x G V and A G DIR, INF(x, A, t) is bimodal in t. It is strictly 
bimodal if and only if 
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(i) ||x|| oo > 6 and 

(ii) -A-lXi > A 2 x 2 > {-AiX! + l)/2 or -A 2 x 2 > A^xi > (-A 2 x 2 + l)/2. 



Proof. A chip at vertex x requires at least time steps to arrive at the origin. 

Hence, INF(x, A, t) = for t < \\^.\\oo- We show that inf(x, A, •) has at most two 
extrema larger than ||x||oo. The discrete derivative of inf(x, A,t) in t is 

p(x,A,t) ■ Ut- l)!) 2 

INF(x,A,t + 2) -iNF(x,A,t) - 



|t+2 ^ t+xr+2 ^ | ^ t-xi+2 ^ | ^ t+x 2 +2 ^ | ^ t-x 2 +2 \ | 

with p(x, A, t) := (4AiXi + AA 2 x 2 )t A + {-Aix\ - A 2 x\ - A x x x x\ - A 2 x 2 x\ - 6A 1 x 1 A 2 x 2 + 
l9A 1 x 1 + 19A 2 x 2 )t 3 + (A 1 x'fA 2 x 2 + A x x x A 2 x\ - AA x x\ - AA 2 x\ - ^A x x x x\ - AA 2 x 2 x\ - 
23A 1 x 1 A 2 x 2 + 30A 1 x 1 + 30A 2 x 2 )t 2 + (Aix\x 2 2 + A 2 x\x\ + AA x x\ A 2 x 2 + AA x xiA 2 x\ - 
AAxxl - 4A 2 xl - AAxXiX?, - AA 2 x 2 xj - 32A 1 x 1 A 2 x 2 + I6A1X1 + lQA 2 x 2 )t - Aix\A 2 xl + 
AAix\A 2 x 2 + AAiX\A 2 x 2 — 16Aix±A 2 x 2 . We observe that the number of extrema of 
inf(x, A, •) is exactly the number of roots of p(x, A, •). Since this a polynomial of 
degree 4 in t, we can use Descartes' Sign Rule and some elementary case distinctions 
to show that p(x, A, •) has at most two roots larger than Hx^. A closer calculation 
reveals that p(x, A, ■) has precisely two roots larger than Hx^ if ||x||oo > 6 and one of 
-AxX! > A 2 x 2 > (-AiXi + l)/2 and -A 2 x 2 > A x x x > (~A 2 x 2 + l)/2 hold. □ 

Lemma 7. For all x e V and A (1) ,A (2) e dir such that A (1) ^ A 2) , 
inf(x, (A^, A^),t) is unimodal in t. 

Proof. The discrete derivative of inf(x, (A^\ A^'),t) is 

inf(x, (A (1) , A (2) ), t + 2) - inf(x, (A (1) , A (2) ), t) 
(p(x,AW,t)+p(x,A 2 ),t))-((t-l)!) 2 

^ t+xi±2 j | ^ t-X!+2 \ | ^ t+x 2 +2 \ | ^ t~x 2 +2 \ | 

with p(x, A, t) as defined in the proof of Lemma EJ As there, the extrema of INF are 
the roots of the quartic function p(x,A^\t) + p(x, A 2 \t). Descartes' Sign Rule now 
shows that p(x, A^\t) + p(x,A^ 2 \t) has at most one root larger than \\x\\oo for all x 
and A 1 ) ^ A 2 ). □ 
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6 Maximal contribution of a vertex 



We now fix a position x and a rotor sequence 7Z to examine maxcon(x). Lem- 
mas [6] and [7J show that ^ Ag ^ inf(x, A, t) is bimodal in t for all finite sequences 
A := (A^, A^, . . .) of rotor directions ordered according to 1Z. Hence, for all A there 
are at most two times at which the monotonicity of ^ Ag ^ inf(x, A, t) changes. A 
time t at which the monotonicity of inf ( x ; A, t) changes for some A is called ex- 

tremal. In case of ambiguities, we define the first such time to be extremal. That is, for 
unimodal ^2 AeA inf(x, A, t), we choose the first time t\ such that ^2 A£A inf(x, A, t) 
is monotone for t < ti and t>t\. Analogously, for strictly bimodal inf ( x ; A, t), 

we choose the first times ti and t 2 such that J2a^a inf ( x > A, t) is monotone for t < t±, 
ti < t < t 2 , and t > t 2 . The set of all extremal times is denoted by ex(x). 

ex(x) can be computed easily. By equation (j4j) it suffices to consider A of length one 
and two. The corresponding extremal times are the (rounded) roots of the polynomials 
p(x, A, t) and p(x, A (1) ,t) +p(x, A( 2 ), t) given in Lemma El The following lemma shows 
that the number of extremal times is very limited. 

Lemma 8. |ex(x)| < 7 

Proof. According to Lemma El there is at most one rotor direction A for which 
inf(x, A, t) is strictly bimodal in t. Hence, the number of extremal times of inf(x, A, t) 
with \A\ = 1 is at most five. For a rotor sequence 1Z = (R^, R^, W- 3 \ R^), equa- 
tion (El) and Lemma [7] show that inf(x, (R (1) , R (2) ), t) = -inf(x, (R (3) , R (4) ), t) and 
inf(x, (R( 2 ),R( 3 )),t) = -inf(x, (R( 4 ),RW),t) are unimodal in t. Therefore, the total 
number of extremal times of inf(x, (A^\ A^), t) with (A^, A^ 2 )) obeying 1Z is at most 
two. □ 

Between two successive times t 1 ,t 2 G ex(x) U {0,T}, ^ Ag ^ inf(x, A, t) is monotone 
in t for all A. Such periods of time [ti, t 2 ] we call a phase. Note that Y2aga inf ( x ; A, t) 
could also be constant in a certain phase. This implies that it is monotonically increas- 
ing as well as monotonically decreasing. To avoid this ambiguity, we use the terms in- 
creasing and decreasing (in contrast to monotonically increasing and decreasing) based 
on the minima and maxima at extremal times ex(x), which are unambiguously defined 
and alternating. We now define precisely when a function ^2 AeA inf(x, A, t) is increas- 
ing or decreasing. Consider the set E of the extremal times of inp( x , A, t) as de- 
fined above. By Lemmas [6] and [7] we know that \E\ e {1, 2}. We call ^2 A&A inf(x, A, t) 
increasing at t if it has a minimum at the maximal t' e E with t' < t or a maximum 
at the minimal t' e E with t' > t. Analogously, we call Y2a&a inf(x, A, t) decreasing 
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at t if has a maximum at the maximal t' £ E with t' < t or a minimum at the minimal 
t' e E with i' > t. 

By abuse of language, let us say that x sends odd chips at time t if /(x, T — t) ^ 
(mod 4). 

Lemma 9. Let C\ be an arbitrary configuration with run-time T > maxEX(x) and let 
CONc 1 (x) be the corresponding contribution o/x. Then there is a configuration C2 with 
the same run-time and C0Nc 1 (x) < CONc 2 (x) that sends odd chips only at extremal 
times, i.e., for the associated f satisfies /(x, T — t) ^ (mod 4) only if t £ ex(x). 

Proof. Let C2 be a configuration with CONc 2 (x) > CON^ (x) and a minimal number of 
non-extremal times at which odd chips are sent from x. We assume this number to be 
greater than zero and show a contradiction. 

The sum of the infs of all chips sent at a certain non-extremal time t is either increasing 
or decreasing in the phase t lies in. 

Let us first assume that it is increasing. Let t' be the minimal t' such that t' £ EX(x) 
or there are odd chips sent at time t' (assume for the moment that such a t' exists). 
Then, sending the considered pile of odd chips at time t' instead of time t decreases the 
number of non-extremal times while not decreasing its contribution. Such a modified 
configuration exists by Theorem [T] and contradicts our assumption on C 2 - Therefore, 
there is no such time t! . This implies that t lies in the last phase and that the odd 
chips sent at time t are the last to be sent at all. By lim^oo inf(x, A, t) = for all A, 
the contribution of the chips sent at time t is negative (since increasing). Hence, not 
sending these chips at all does not decrease CONc 2 (x), but the number of non-extremal 
times. 

The same line of argument holds if the sum of the infs is decreasing instead of increas- 
ing. In this case we use that inf(x, A, t) = for all t < ||x||oo. □ 

Lemma [9] immediately gives the following corollary. 

Corollary 10. There is an initial configuration and a time T such that CON(x) = 
maxcon(x). The configuration can be chosen such that /(x, T — t) ^ (mod 4) only 
if t £ ex(x). T can be chosen arbitrarily as long as T > maxEx(x). 

Lemma El and Corollary \IU\ already give a simple, but costly approach to calculate 
MAXCON(x): There are four different initial rotor directions for x and at each (of 
the at most seven) extremal time we can either send 0, 1, 2, or 3 odd chips. As 
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all subsequent rotor directions are chosen according to TZ, there are only a constant 
4 • 4 7 = 65536 configurations to consider. The maximum of the respective con(x) will 
be MAXCON(x) by Corollary [TOl 

Fortunately, we can also find the worst-case configuration directly. A block of a phase 
[*i,t 2 ] is a 4-tupel (A (1) , A (2) , A (3) , A (4) ) G dir 4 of rotor directions in the order of TZ 
such that Yli=i lNF ( x ; AW, t) is increasing in t in this phase for all k G {1,2,3}. By 
equation ([3]), this is equivalent to Ylt=k INF ( X > A^, t) being decreasing in t within the 
phase for all k G {2, 3, 4}. 

Lemma 11. Each phase has a unique block. This is determined by the monotonicities 
of inf(x, A,t) with \A\ G {1,2}. 



Proof. Consider a fixed phase. We want to show that for all valid combinations of 
monotonicities of inf(x, A, t) with \A\ G {1,2} within this phase, there is exactly one 
permutation (A (1) , A (2) , A (3) , A (4) ) of dir obeying TZ such that (A (1) , A (2) , A (3) , A (4) ) 
forms a block. 

To describe the type of monotonicity of inf(x, A, t) within the phase, we use a function 
r with t(A) :=— > if inf(x, A,t) is increasing and t(A) :=<— if it is decreasing. This 
notation should indicate the direction in which the respective inf(x, A, t) is increasing. 
As a short form we also use r (AW, A( 2 ), AW, AW) := (t(A^),t(A^),t(A^),t(A^)). 

By equation ([3]), we know that there is at least one A of type — If there is exactly 
one direction A of type then the unique permutation (A^, A^ 2 \ A^ 3 \ A^) of dir 
obeying TZ such that r(A^\ A> 3 \ A^ 4 ') = (—►,«—,•<—,•<—) is the uniquely defined 
block. If there are three rotor directions A of type — >, the block is analogously uniquely 
defined by r(A^, A®, A^ 3 \A^) = (->, <-). 

It remains to examine the case of exactly two rotor directions of type — If these two 
directions are consecutive in TZ, r(A^\ A^ 2 ', A- 3 ', A^) = (— >, — >, <— , ) again defines 
the unique block. Otherwise, rotor directions of type — > and <— are alternating in the 
rotor sequence and (— >, <— , — >, <— ) is the only type possible for a block. This allows two 
blocks (A (1) , A (2) , A (3) , A (4) ) and (A (3) , A (4) , A (1) , A (2) ). The choice between these two 
is uniquely fixed by the monotonicity of inf(x, (A^\ A^),t). Therefore, in all cases 
there is exactly one unique block. □ 



We now use Lemma [TT] to define a particular configuration, which we call block con- 
figuration. By Theorem [T], we may specify a configuration sufficiently well by fixing 
the number of odd chips at all times and locations. In a block configuration B, a ver- 
tex x sends odd chips only at extremal times t G ex(x). Let (M 1 ', A^ 2 ', M 3 \ A^) and 
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Figure 1: iNF((g), A,t) for A € {/",\, The circles indicate the extrema. 



(A^), A'- 2 '', A^ 3 \ A^) denote the blocks in the phases ending and starting at t. Then x 
sends k chips at time t in directions (A^ 1 ), . . . , A^), where k is such that < k < 3 
and (. . . , A^, A^), . . . , A^ k \ A^\ . . .) obeys 1Z. This uniquely defines when and in 
which directions odd chips are sent. Note that we used the blocks only as a technical 
tool. There are not necessarily chips sent corresponding to A^ , AS 2 \ M 3 ' , A^ and 
AS 1 ', M 2 \ A® , A( 4 ). By Theorem (TJ there are configurations B as just defined and for 
all x all of them have the same contribution cONe(x). 

Example. We now derive the block configuration of the position x = with the clock- 
wise rotor sequence 1Z = (/*, \, y , \). By calculating the roots of the polynomials 
p(x, A, t) and p(x, A^ , t) + p(x, A^ 2 ) , t) given in Lemma El it is easy to verify that 

• INF(x, /*, t) is unimodal with minimum at t = 27. 

• inf(x, \, t) is bimodal with minimum at t — 9 and maximum at t = 35, 

• inf(x, y , t) is unimodal with maximum at t = 25, 

• INF(x, \, t) is unimodal with minimum at t = 23, 

• inf(x, (/*, \), t) and inf(x, (\, /*), t) are unimodal with minimum at t = 27. 

• INF(x, (\, S), t) and INF(x, (,/, \), t) are unimodal with maximum at t = 27. 

Hence, the extremal points are ex(x) = {9,23,25,27,35}. Figure [j] depicts the plots 
of inf(x, A, t). The modes of inf(x, A, t) listed above uniquely determine the blocks 
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of each phase. The following table lists rotor directions and type of the block of each 
phase. 



Boundaries of the phase Block of the phase 

lower upper Rotor directions Type 

9 

1 9 23 

2 23 25 

3 25 27 

4 27 35 

5 35 T 



This yields the following (maximal as we will see shortly) contribution at x = : 

con(x) = inf(x, /, 9) + inf(x, \, 9) + inf(x, /*, 9)+ 
inf(x,\,27) + inf(x,/,27) 

20,506,216,364,597 ^ q nf)9977 
~~ 9,007,199,254,740,992 ~ U.UUZZ (( . 

Note that just sending a single chip in the worst direction y at its worst time t = 25 
gives a smaller contribution of inf(x, 25) ~ 0.001985. Also, sending two chips in 
directions \ and / at time 27 = argmax t inf(x, (\, ,/),t) gives inf(x, (\, >/), 27) ~ 
0.002261. Hence we do profit from sending a chip in the "wrong" direction /* at time 
9. 

The values of CON ( ( 9 ) ) for other rotor sequences are shown in the following table. 



Times and directions of odd ,, 5 
Rotor sequence chips in a bbck configuration CON ( ( ft 



(/,\,/,\) 9: 


/\/, 27 :\/ 


0.002277. . . 


(/,\,/,\) 23:/^ 


\A 27 :\/, 35 :\ 


0.002309. . . 




, 23:\//\ 27:\\ t 


/ 0.002302... 




25 35 :\ 


0.002230. . . 




27:\\/ 


0.002083. . . 




25 :/ 


0.001985... 



Lemma 12. A block configuration yields a contribution o/MAXCON(x). 

Proof. Consider a configuration C with contribution cONe(x) = maxcon(x). By pre- 
vious considerations, we can further assume the following. 
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(1) C only sends odd chips at times t G ex(x) (cf. Corollary [TUj) . 

(2) C sends at least seven chips at each time t G ex(x) (cf. equation ([3])). 

(3) Let ti,t 2 G ex(x) such that [ti,t 2 ] is a phase and let k G {1,2,3}. Let 
A^, . . . , A( fc ) be the directions the last k chips are sent from vertex x at time t\. 
If Yli=i INF ( X ; AW, t) is increasing (cf. definition on page [12]) in [ti,^], then it is 
not constant. This is a feasible assumption on C, since otherwise we could send 
these k chips at time t 2 without changing cONc(x). 

(4) Analogously, let ti,t 2 G ex(x) such that [ti,t 2 ] is a phase and let j G {1,2,3}. 
Let A 1 ), . . . , A^ be the direct ions the first j chips are sent from vertex x at time 
t 2 . If Yli=i lNF ( x ; AW, t) is decreasing in [£1,^2], then it is not constant. 

Let B be a block configuration. Aiming at a contradiction, we assume CONc(x) > 
CONe(x). Since by Assumption (CD) and the definition of B both configurations send 
odd chips only at times in EX(x), there is a time t G EX(x) at which the chips of C 
contribute more than the chips of B. 

We now closely examine the chips sent from x at time t by both configurations. We 
know that B sends a uniquely determined number £ G {0, . . . , 3} of odd chips at time t 
in some directions A^, . . . , A^. By the above Assumption ([2]), C also sends a sequence 
of chips in directions A^, . . . , A™. Let j and k denote the number of chips sent by C 
at time t before and after these I chips, respectively. By ignoring, possible piles of four 
chips, we may assume j, k < 3. 

Assume that k > 1 . Then the sum of the infs of the last k chips C sends at time t is 
increasing by the definition of a block. Assume first that t is not the last extremal time, 
that is, there is some t 2 G ex(x) such that [t, t 2 ] form a phase. Then by Assumption ([H]) 
above, the sum of the infs of the last k chips is strictly increasing in [t,t 2 ]. Hence, 
a configuration which sends these chips instead at t 2 has a larger contribution, in 
contradiction to the maximality of C. Now let t be the last extremal time. From 
lim^oo inf(x, A, t) = for all A and the fact that the sum of the infs of the last 
k chips is increasing, we see that it is not positive. Hence the last k chips do not 
contribute positively to CON c (x). 

Analoguously, assume that j > 1. Assume first that t is not the first extremal time, that 
is, [ti,t] form a phase for some ti G ex(x). By Assumption (jl]), the first j chips C sends 
at time t have a strictly monotonically decreasing sum of infs. Hence sending them 
at time t± instead of t gives a larger contribution, again contradicting the maximality 
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of C. If t is the first extremal time of x, then inf(x, A, t) = for all A and t < \\^\\oo 
shows, similarly as above, that the contribution of the first j chips is not positive. 

We conclude that the first j and last k chips sent from x and time t in C, if they are 
present, do not contribute positively to the contribution of x. This contradicts our 
assumption CON c (x) > CON B (x). □ 

With the help of a computer, we can now calculate maxcon(x) for all x. Using about 
two months on a Xeon 3 GHz CPU, we computed the maximal contribution of all 
vertices in [—800, 800] 2 . If we have the same rotor sequence for all vertices then 



E 

:||oo<800 



. 7.832... for a circular rotor sequence 
maxcon(x) = \ (5) 
7.286... for a non-circular rotor sequence. 



On the other hand, if we allow a different rotor sequence for each vertex, and further 
assume that each vertex has a rotor sequence leading to the maximal contribution, 
then we get 

maxcon(x) = 7.873... 

Hxlloo^SOO 

Since all maxcons are non-negative, above values are a lower bound for 
Sxgv max COn(x), and hence for the single vertex discrepancy by Theorem HI 

Remark. Lemma [8] shows that the number of extremal times of a vertex is at most 
seven. However, a block configuration does not send odd chips at at all extremal times. 
Let ex(x) denote the set of extremal times at which odd chips are sent by the block 
configuration. Figure H] at the end of the paper shows that there are vertices x such 
that |ex(x)| > 3. We now sketch a proof that |ex(x)| < 4 for all x. 

Note that |ex(x)| only depends on the relative order of the extremal points and the 
initial monotonicity (i.e., increasing or decreasing) of INF(x, A, t) for |*4| < 2. We use 
the following two properties of INF (derived from equation Q): 

• In each phase there is at least one A e dir such that inf(x, A, t) is increasing 
(or decreasing). 

• If inf(x, A (1) ,t) and inf(x, A (2) , t) are both increasing or decreasing in a phase, 
so is inf(x, (A 1 ), A 2 )),t). 

For a vertex x with only unimodal inf(x, A,t), there are 6! = 720 permutations of 
the extrema of inf(x, A, t) and inf(x, (A^, A^),t) and 2 6 = 64 initial monotonicities 
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(using equation (j3D). A simple check by a computer shows that for only 384 of these 
46080 cases both properties from above are satisfied. For all of them, |ex(x)| < 3 holds. 
For vertices x with inf(x, A, t) strictly bimodal for an A £ dir, there are 7\/2\ = 2520 
permutations of the extrema and 2 6 = 64 initial monotonicities. Here, all 408 cases 
which satisfy both properties only achieve |ex(x)| < 4. This proves |ex(x)| < 4 for 
all x. Figure H] shows |ex(x)| < 3 for all < 10. We could also verify this for 

|| x l|oo < 800. Therefore, we actually expect |ex(x)| < 3 to hold for all x. To bridge 
this gap, stronger properties of INF seem necessary. 



7 Tail Estimates 

In the previous section, we have calculated the values of Xl||x|| oo <800 maxcon(x) de- 
pending on the rotor sequence. To show that these are good approximations for the 
maximal single vertex discrepancy, we need to find an upper bound on 

E := MAXCON(x). 

||x|| oo >800 

In this section, we will prove E < 0.16. 

We now fix an arbitrary initial configuration and a time T. A simple calculation based 
on the definitions of INF and CON gives for all x, A, and t 

inf(x, A, t) = ((A lXl ■ A 2 x 2 )t~ 2 - (Aixi + A 2 x 2 )t~ 1 ) if(x, t), 

Si(x)<T 

with Sj(x) as defined in Section 0] and A^ := next 1 (arr(x, 0)). Note that, inde- 
pendent of the chosen rotor sequence, each of the sequences (A^\x))i> , {A\ ' (x))j>o, 
and (Af(x)Af(x)) i> is alternating or alternating in groups of two. To bound the 
alternating sums in equation ©, we use the following fact, which is an elementary 
extension of Lemma 4 in . 

Lemma 13. Let f'X — > M. be non-negative and unimodal with X C R. Let 
A (0) , . . . , A {n) e {-1,+1} andt ,...,t n £ X such that t < ... < t n . If A® is al- 
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ternating or alternating in groups of two, then 



n 



$> (0 /(*i) < 2™:/(x). 



It remains to show that H(x,t)/t and i?(x, t)/t 2 are indeed unimodal. Note that 
inf(x, A, t) itself is not always unimodal as shown in Lemma EJ 

Lemma 14. For all x £ V, H(?t,t)/t and H(x,t)/t 2 are unimodal in t with global 
maxima at t max (x) and t' m3X (x), respectively. For the maxima {x\ + x 2 )/A — 2 < 
tmax(x) < (x\ + x 2 2 )/A + 1 and (x\ + z|)/6 - 1 < C ax ( x ) < ( x i + X l)/ G + 2 ^ oWs - 

Proof. By symmetry, let us assume X\ < Xi- By definition, H(x.,t)/t = for t < Xi- 
We show that H(jc,t)/t has only one maximum in t £ [x 2) oo). We compute 



withp(t) := At 3 — {x\ + x\ + h)t 2 + 2t + x\x\. By Descartes' Sign Rule (cf. Theorem^)), 
p(t) has at most one real root larger than x 2 - Since 



we see that H(x,t)/t has a unique extremum, which is a maximum, in \{x\ + x 2 )/A — 
2, (x^ + x 2 )/A + 1]. This proves the lemma for H(x,t)/t. The analogous proof for 



x,t-2) #(x,t) 
i - 2 i 



4"V(*)(*-3)! 2 (t-2) 




i7(x, t)/t 2 is omitted. 



□ 



By equation (jSJ), Lemmas [T31 and [T41 we obtain 



E < AE X + 2£ 2 



with 




Using Lemma [Hand #(x,t) < (2 _ *( t * 2 )) < 1/f, we now derive upper bounds for 
H(x,t)/t and H(x,t)/t 2 for (jx^ > 88: 
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#(X, tmax(x)) 



^max (xj 



< 



^max(x) 



#( X , Cax(x)) 



< 



^max(x) 



< 



< 



16 



< 



17 



(x 2 + :r 2 -8) 2 - {x\ + xl) 2 ' 
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(x 2 + X 2 . — 6) ; 



< 
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For the calculations in the remainder of this section we need the following estimates. 
All of them can be derived by bounding the infinite sums with integrals. 



Y 1 



< 



1 



x>y 



x k (k — l)y k 



— for all y > and all constants k > 1. 



oo 1 7 

51 — 2T2 — o - ^ forallxi>l. 



V 1 

/7v2 4- ■ " 



(a 2 + 1/ 2 )|/ 



E 



ln(a 2 + /3 2 ) -21n(/3) 
2^ ' 

1 

< 



y>a, 
y=c(mod 2) 



(?/ 2 + 7 2 ) 2 - 4(a 2 + 7 2 ) 



E 



1 ^ (?r - 2arctan(^))(a 2 + 7 2 ) - 2«7 



y>a, 
y=c(mod 2) 



(?/ 2 + 7 2 ) 5 



(a 2 + 7 2 )7 3 



E 



2arctan(|) hi(a 2 + /3 2 ) -21n(/3) , tt - 2 arctan(^) 



ST 



+ 







With this, we can now bound E2 easily: 

217xiX 2 



e 2 < y: 



< 



||x|| oo >800 
800 

E E 



{x\ + x 2 .) 3 



£ E 

||x|| cx) >800 



217 



2(x? + x 2 ) 2 



434 



xi=i x 2 >soo ( Xl + x ' 2) 

800 



2\2 



^ x ^ 434 

1^ 2^ ( x 2 + x 2)2 



Ex - 434 x - 

J 9 „ „ „ \ 



a'!>800 x 2 >0 

3038 



< 



xi=lx 2 >800 2 m>800 

434 1519 



3 • 800 2 3 • 800 2 



< 0.0011. 



(7) 
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Achieving a good bound for E\ is significantly harder. We divide E\ in three subsums: 

see equation J9} see equation IllOjl 
^800 A * ' 800 " 

^ < 4 XiH(x, t max (x)) ^ ^ xiH(x, t max (x)) 

xi=l xj,=801, *max(x) t m ax(x) 

I2=a;i (mod 2) X2=xi(mod 2) 



oo 



+ 4 £ £ Xlg(M - (x)) < 0.038. 

B1 =801 xa=801, ^maxlXj 
X2=a;i(mod 2) 



see equation (1111) 

Now we bound these sums separately as follows. 

800 oo IJ{ + [ x\ 800 

4 ^ ^ ^(x,t (x)) < 6g ^ xi 



tmax ( X ) xi=l x 2 >800, ( X l + X 2) 2 

/T2=a;i (mod 2) £2=xi(mod2) 

17 |™ (800 2 + x*) (tt - 2 arctan(800/xi)) - 1600 x l 

< — V ^ — ^ ' J < 0.0046. 9 

~ 2 ^ x? + 800 2 )x\ V ' 

oo 800 rr/ . / \ \ 800 

4 ^ £ Xl max( - X ^ < 68 E £ 



^max(x) „„. (^i ~~\~ Xn) 

a;i=801 X2=0, v ' £2=0 xi>800, v 1 

X2=xi (mod 2) xi=X2(mod2) 



< 17 < 0.0167 (10) 



— ( n 

22=0 2 



4 g g x x i7(x,t (x)) ^ 6g ^ ^ ^ 



^max( x ) xi>800 x 2 >800, ~^ X ^ 



< 



i2=ii(mod 2) X2=xi(mod2) 

17 (tt -2arctan(^))(800 2 + x 2 ) - 1600xi 

T ^ (800 2 + x?)x? 

xi>800 v '1/1 



17 f\n(2 ■ 800 2 ) - 21n(800) , tt/2 ln(2 ■ 801 2 ) - 2 ln(801) 



2 V 800 800 801 

< 0.0167. (11) 

Putting this together, we obtain 

E < 4- 0.038 + 2 ■ 0.0011 < 0.16. (12) 

This upper bound on E is not tight. However, it suffices to prove that the bounds for 
the single vertex discrepancy calculated in Section [6] do depend on the rotor sequence. 
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Theorem H] and Equations (jSJ) and (fT2*|) yield the following theorem. 



Theorem 15. TTie maximal single vertex discrepancy between Propp machine and 
linear machine is a constant ci, which depends on the allowed rotor sequences: 

• If all vertices have the same circular rotor sequence, 7.832 < C2 < 7.985. 

• If all vertices have the same non-circular rotor sequence, 7.286 < C2 < 7.439. 

• // all vertices may have different rotor sequences, and we assume that each vertex 
has a rotor sequence leading to a maximal contribution, then 7.873 < c% < 8.026. 



8 Aggregating Model 

Besides the small single vertex discrepancies examined in the previous sections, Propp 
machine and random walk bear striking similarities also in other respects. The histori- 
cally first research started by Jim Propp regarded an aggregating model called Internal 
Diffusion-Limited Aggregation (IDLA) [sj]. In physics this is a well-established model 
to describe condensation around a source. 

The process starts with an empty grid. In each round, a particle is inserted at the origin 
and does a (quasi)-random walk until it occupies the first empty site it reaches. For 
the random walk, it is well known that the shape of the occupied locations converges 
to a Euclidean ball jfl in the following sense. Let n be the number of particles and 
let A(n) denote the difference of the radius of the largest inscribed and the smallest 
circumscribed circle of an aggregation with n chips. It has been shown by Lawler 
\^ that the fluctuations around the limiting shape are bounded by @{n l / G ) with high 



probability. Moore and Machta lOj observed experimentally that these error terms 



were even smaller, namely poly-logarithmic. 

The analogous model in which the particles do a rotor-router walk instead of a random 
walk is much less understood. Levine 0] showed that after n particles have been 
added, the Propp aggregation contains a disc of radius Vt{n 1 ^). Levine and Peres [si, 0] 
proved that the shape of occupied locations converges to a Euclidean ball, however, in 
a weaker sense than before. They showed that the Lebesgue measure of the symmetric 
difference between the Propp aggregation and an appropriately scaled Euclidean ball 
centered at the origin is (^(n 1 / 3 ). Surprisingly, experimental results indicate much 
stronger bounds. Kleber 0] computed that for counter-clockwise permutations of the 
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(a) Counterclockwise rotor sequence (f , *—, j, — >). 



I: . ^ s 




(b) Non-circular rotor sequence (t, *—, — >, !)■ 



igure 2: Propp aggregations with one million particles. All rotors initially point to the 
left. The final rotor directions up, left, down, and right are denoted by the colors red, 
yellow, green, and blue, respectively. 
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rotor directions A(3 • 10 6 ) ~ 1.611 if all rotors initially point to the left. An apparent 
conjecture is that there is a constant S such that A(n) < 5 for all n. 

We reran these experiments with different rotor sequences. The aggregations for one 
million particles are shown in Figure El Both aggregations do not only differ in the 
color patterns, but also in the precise value of A(n). If all rotors are initially set to the 
left, we obtained the following values for A(n). 



Rotor sequence (<-, T,->,|) (<-,T>Jo-0 



average A(n) 


for 2 • 


10 6 < n < 3 • 10 6 


1.600 


0.996 


1.810 


maximal A(n 


) for n 


< 3 • 10 6 


1.741 


1.218 


1.967 



It is noteworthy that the respective A(n) of both non-circular rotor sequences («— , f , I 
, — >) and («—,—>•, t> I) differ considerably. 

Additionally, we also examined A(n) for random initial rotor directions. This leads to 
slightly larger A-values. The following table shows averages and standard deviations 
of 100 aggregations with random initial directions of the rotors. 



Rotor sequence (<-, t, I) (<-»t>-U-Q ~S t,l) 

average A(n) for 2 • 10 6 < n < 3 • 10 6 1.920 ± 0.004 1.782 ±0.003 1.781 ±0.003 
maximal A(n) for n < 3 • 10 6 2.541 ± 0.051 2.351 ± 0.053 2.364 ± 0.067 



As one might have expected, for random initial rotor directions the two non-circular 
rotor sequences (columns one and three) are statistically not distinguishable. 

The results above again show that different rotor sequences do make a difference. The 
main open problem, however, remains to show the conjectured constant upper bound 
for A(n). 



9 Conclusion 



One way of comparing the Propp machine with a random walk is in terms of the 
maximal discrepancy that can occur on a single vertex. It has been shown by Cooper 
and Spencer [l| that for the underlying graph being an infinite grid Z d , this single 
vertex discrepancy can be bounded by a constant independent of the particular 
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initial configuration. For d = 1, this constant has been estimated as C\ ~ 2.29 in [2|. 
Also, the initial configurations leading to a high discrepancy have been described. For 
d > 2, no such results were known. 

In this paper, we regarded the case d = 2. We chose the case d = 2 out of two 
considerations. On the one hand, from dimension two on, there is more than one 
rotor sequence available, which raises the question if different rotors sequences make 
a difference. One the other hand, we restrict ourselves to d — 2, because for larger 
d a nice expression for the probability -ff(x, t) that a chip from vertex x arrives at 
the origin after t random steps is missing. This probably makes it very hard to find 
sufficiently sharp estimates for the single vertex discrepancies. 

We were able to give relatively tight estimates for the constants c 2 taking into account 
different rotor sequences and obtain several interesting facts about the worst-case initial 
configurations. The maximal single vertex discrepancy c 2 satisfies the following. If all 
vertices have the same circular rotor sequence, 7.832 < c 2 < 7.985. If all vertices 
have the same non-circular rotor sequence, 7.286 < c 2 < 7.439. If all vertices may 
have different rotor sequences, and we assume that each vertex has a rotor sequence 
leading to a maximal contribution, then 7.873 < c 2 < 8.026. In particular, we see 
that non-circular rotor sequences seem to produce smaller discrepancies than circular 
one. The gaps between upper and lower bounds stem from the fact that we used a 
computer to calculate the precise maximal contribution CON(x — y) of vertex x on 
the discrepancy at y. Hence the lower bounds are the maximal discrepancies obtained 
from initial configurations such that all vertices x with ||x — y||oo > 800 at all times 
contain numbers of chips only that are divisible by 4. 

We also learned that the initial configurations leading to such discrepancies are more 
complicated than in the one-dimensional case. Recall from that in the one- 
dimensional case in a worst-case setting each position needs to have an odd number of 
chips only once. If we aim at a surplus of chips in the Propp model, these odd chips 
were always sent towards the position under consideration, otherwise away from it. 

In the two-dimensional case, things are more complicated. Here it can be necessary 
that a position holds a number of chips not divisible by 4 up to three times. Also, the 
number of "odd" chips (those which cannot be put into piles of four) can be as high as 
nine. In consequence, it can make sense to send odd chips in the wrong direction (e.g., 
away from the position where we aim at a surplus of chips). An example showing this 
was analyzed in Section [6l The reason for such behavior seems to be that the influences 
inf(x, A, t) of odd chips sent from x in direction A at time t are not unimodal functions 
in t anymore (as in the one-dimensional case). 
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In Figures [H] to El more information about the behavior of different positions in a 
worst-case setting (aiming at a surplus of chips at the origin) is collected. 

We also briefly regarded the IDLA aggregation model. We saw that the surprisingly 
strong convergence to a Eucledian ball observed in earlier research also holds for non- 
circular rotor sequences and non-regular initial rotor settings. However, the suspected 
constant again seems to depend on the rotor sequences, and again, the circular ones 
seem to behave slightly worse than the non-circular ones. 
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(b) Non-circular rotor sequence (/*, \, \, /). 

Figure 3: At each position, the directions in which odd chips are sent in a block configura- 
tion are shown. If there are multiple chips sent from one position, the arrows are scaled 
according to their relative contribution. Grey arrows represent negative contributions. For 
position ( g ) , the second and third row of the second table on page [16] contains more details 
on the times when these chips are sent. on 



1 1 10i 2 3 



2 2 2 1 1 
3 2 2 1 1 

3 2 2 1 1 
3 3 2 2 1 

3 3 2 2 15 
3 3 3 2 1 

3 3 3 1 1 
2 2 3 2 1 

11111 



~~ I X — X — 1 X — TT 

- 10 1 1 "I 3 1 I 
11111 

112 11 
112 2 2 

2 2 2 3 15 
2 2 2 3 3 

2 2 3 3 1 
2 2 3 3 2 

2 3 3 3 1 



1 3 3 3 2 

2 3 3 2 2 
1 3 3 2 2 

3 3 2 2 2 
1 3 2 2 2 

2 2 2 1 1 
112 11 

11111 
11111 
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11x11 

1 2 3 2 2 
113 3 3 

1 2 3 3 3 
1 2 2 3 3 

1 2 2 3 3 
112 2 3 

112 2 3 
112 2 2 
3 3 3 2-1 0^ 112 2 2 

(a) Counterclockwise rotor sequence (/*, \, S , \). 
2222 10i 12332 



1 2 2 2 1 
112 2 2 

112 2 1 
1112 2 

1 1 1 2 15 
11112 

11111 
11111 

11111 



~I 1 — 1 — 1 X — TT 

- 10 1 1 "I 3 1 I 
11111 
11111 

11112 
1 1 1 2 15 

1112 2 
112 2 1 

112 2 2 
1 2 2 2 1 



1 2 3 3 2 

2 2 3 2 2 
1 2 3 2 2 

2 3 2 2 2 
1 3 2 2 2 

2 2 2 2 2 
1 2 2 2 2 

1 2 2 2 2 
112 2 2 
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1 2 2 2 2-10} 1 2 

(b) Non-circular rotor sequence 



112 2 2 

1 2 2 2 2 
1 2 2 2 2 

2 2 2 2 2 
1 3 2 2 2 

2 3 2 2 2 
1 2 3 2 2 

2 2 3 2 2 
1 2 3 3 2 

3 3 2 



5 \7 i 



Figure 4: Number of times at which a block configuration sends odd chips. Note that the 
remark on page 1181 only proved |ex(x)| < 4, though there is no position x observable with 
|ex(x)[ > 3. 
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5 5 5 2 2 102 

5 5 5 2 2 

6 5 5 2 2 

6 5 5 2 2 
6 6 5 5 2 

6 6 5 5 25 
6 6 6 5 1 

6 6 6 1 2 
2 6 6 2 1 

2 2 2 2 1 



~z — — z — — z — i — z — — z — -Q- 
_10 2 2 ^2 1 

2 2 2 1 1 

2 2 5 1 2 

2 2 5 5 2 

5 5 5 6 25 

5 5 5 6 6 

5 5 6 6 2 

5 5 6 6 6 

5 6 6 6 2 

5 6 6 6 2-102 



2 6 6 6 5 
2 6 6 6 5 

6 6 6 5 5 
2 6 6 5 5 

6 6 5 5 5 
2 6 5 5 5 

2 5 5 2 2 
2 15 2 2 

112 2 2 
1 2 2 2 2 
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1 2 2 2 2 
1 2 6 6 2 

2 16 6 6 

1 5 6 6 6 
2 5 5 6 6 

2 5 5 6 6 
2 2 5 5 6 

2 2 5 5 6 
2 2 5 5 5 
2 2 5 5 5 



(a) Counterclockwise rotor sequence (/", \, S , \). 



1 2 2 2 2 102 

1 2 2 2 2 
112 2 2 

112 2 2 
1112 2 

1 1 1 2 25 
11112 

11112 
11111 

11111 



x x . x x 

u i i r i i 



11111 

11112 
11112 

1 1 1 2 25 
1112 2 

112 2 2 
112 2 2 

1 2 2 2 2 
1 2 2 2 2-102 



2 6 9 9 5 
2 6 9 9 5 

6 6 9 5 5 
2 6 9 5 5 

6 9 5 5 5 
2 9 5 5 5 

5 5 5 5 5 
2 5 5 5 5 

1 5 5 5 5 
115 5 5 



1 1 I 5 5™ 

1 5 5 5 5 

2 5 5 5 5 

5 5 5 5 5 
2 9 5 5 5 

6 9 5 5 5 
2 6 9 5 5 

6 6 9 5 5 
2 6 9 9 5 

2 6 9 9 5 



(b) Non-circular rotor sequence (/*, \, \, /). 

Figure 5: Number of odd chips sent in a block configuration. The largest number of 
odd chips occurs only in |(b)[ not in (a) 
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0.8 1.1 1.5 2.1 2.6KJ-8 2.6 2.1 1.6 1.1 0.8 

1.1 1.6 2.3 3.1 3.7 3.7 3.1 2.3 1.6 1.1 
1.1 1.6 2.4 3.6 4.8 5^5 4.9 3.6 2.4 1.6 1.1 

1.6 2.5 3.9 5.9 7.9 7.9 6.1 3.9 2.5 1.6 
1.6 2.4 4.1 7.0 11 13 11 7.0 4.1 2.4 1.5 

2.3 3.9 7.4 14 215" 21 14 7.4 3.9 2.3 
2.1 3.6 7.0 16 31 47 34 16 7.0 3.6 2.1 

3.1 6.1 14 47 94 94 47 14 5.9 3.1 
2.6 4.9 11 34 188 2!i0 188 31 11 4.8 2.6 

3.7 7.9 21 94 750 750 94 21 7.9 3.7 

2j8 5:5 13 i 47 250 q 250 ■ 47 i 13 5:5 2j8 

~^.7 7.9 2^ 94 750 750 94 A 7.9 3.7^ 
2.6 4.8 11 31 188 2!i0 188 34 11 4.9 2.6 

3.1 5.9 14 47 94 94 47 14 6.1 3.1 
2.1 3.6 7.0 16 34 47 31 16 7.0 3.6 2.1 
2.3 3.9 7.4 14 215- 21 14 7.4 3.9 2.3 

1.5 2.4 4.1 7.0 11 13 11 7.0 4.1 2.4 1.6 
1.6 2.5 3.9 6.1 7.9 7.9 5.9 3.9 2.5 1.6 

1.1 1.6 2.4 3.6 4.9 5:5 4.8 3.6 2.4 1.6 1.1 

1.1 1.6 2.3 3.1 3.7 3.7 3.1 2.3 1.6 1.1 
0.8 1.1 1.6 2.1 2. 6-1 02- 8 2.6 2.1 1.5 1.1 0.8 

(a) Counterclockwise rotor sequence (/*, \, S , \). 

0.8 1.1 1.5 2.1 2 . 6 1 Od 8 2.6 2.1 1.5 1.1 0.8 

1.1 1.6 2.3 3.1 3.7 3.7 3.1 2.3 1.6 1.1 
1.0 1.6 2.4 3.6 4.9 5:5 4.9 3.6 2.5 1.6 1.1 

1.5 2.4 3.9 6.0 7.9 7.9 6.0 4.0 2.5 1.6 
1.4 2.3 4.0 6.9 11 13 11 7.2 4.2 2.4 1.4 
2.0 3.7 7.2 14 215- 21 14 7.8 3.9 2.1 

1.6 3.0 6.5 16 34 47 33 17 7.0 3.2 1.7 
2.3 4.9 13 47 94 94 52 14 5.1 2.4 

1.6 3.3 8.3 31 188 2!i0 188 33 8.7 3.5 1.7 

2.2 4.8 14 78 750 750 78 15 4.9 2.2 
li4 . 2:7 . G:5 23 ■ 125 . n 125 23 G:5 2:7 1:4 

-10 -5 u 5 10 

2.2 4.8 14 78 750 750 78 f5 4.9 2.2 
1.6 3.3 8.3 31 188 2:i0 188 33 8.7 3.5 1.7 

2.3 4.9 13 47 94 94 52 14 5.1 2.4 
1.6 3.0 6.5 16 34 47 33 17 7.0 3.2 1.7 

2.0 3.7 7.2 14 215- 21 14 7.8 3.9 2.1 
1.4 2.3 4.0 6.9 11 13 11 7.2 4.2 2.4 1.4 

1.5 2.4 3.9 6.0 7.9 7.9 6.0 4.0 2.5 1.6 
1.0 1.6 2.4 3.6 4.9 5:5 4.9 3.6 2.5 1.6 1.1 

1.1 1.6 2.3 3.1 3.7 3.7 3.1 2.3 1.6 1.1 
0.8 1.1 1.5 2.1 2. 6-1 02- 8 2.6 2.1 1.5 1.1 0.8 

(b) Non-circular rotor sequence (/*, \, \, /). 

Figure 6: Contributions MAXCON(x) times 1000. Note the symmetries discussed in Section[6] 
and the quick descent of MAXCOn(x) for increasing 1 1 x| 1 2 - 
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